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It is now widely accepted that the cuprate superconductors are characterized by the same long- 
range order as that present in the Bardeen-Cooper-Schrieffer (BCS) theory: that associated with 
the condensation of Cooper pairs. The author argues that many physical properties of the cuprates 
require interplay with additional order parameters associated with a proximate Mott insulator. 
A classification of Mott insulators in two dimensions is proposed. Experimental evidence so far 
shows that the class appropriate to the cuprates has collinear spin correlations, bond order, and 
confinement of neutral, spin 5 = 1/2 excitations. Proximity to second-order quantum phase tran- 
sitions associated with these orders, and with the pairing order of BCS, has led to systematic 
predictions for many physical properties. In this context the author reviews the results of re- 
cent neutron scattering, fluxoid detection, nuclear magnetic resonance, and scanning tunnelling 
microscopy experiments. 
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The discovery of high temperature superconductiv- 
ity in the cup rate series of compounds by Bednorz and 
Miillcr ( 1986| ) has strongly influenced the development 
of condensed matter physics. It stimulated a great deal 
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of experimental work on the synthesis and characteriza- 
tion of a variety of related intermetallic compounds. It 
also reinvigorated theoretical study of electronic systems 
with strong correlations. Technological applications of 
these materials have also appeared, and could become 
more widespread. 

Prior to this discovery, it was widely assumed that 
all known superconductors, or superfluids of neutral 
fermions such as 3 He, were described b y the theory o f 



Bardeen, Cooper and Schrieffer (BCS) ( Bardeen et al 



1957). Certainly, the quantitative successes of BCS the- 



ory in describing an impressive range of phenomena in 
the lower temperature superconductors make it one of 
the most successful physical theories ever proposed. Soon 
after the discovery of the high temperature superconduc- 
tors, it became clear that many of their properties, and 
especially those at temperatures (T) above the super- 
conducting critical temperature (T c ), could not be quan- 
titatively described by the BCS theory. Overcoming this 
failure has been an important motivation for theoretical 
work in the past decade. 

One of the purposes of this article is to present an 
updated assessment of the applicability of the BCS the- 
ory to the cuprate superconductors. We will restrict our 
attention to physics at very low temperature associated 
with the nature of the ground state and its elementary 
excitations. This will allow us to focus on sharp, qualita- 
tive distinctions. In particular, we will avoid the regime 
of temperatures above T c , where it is at least possible 
that any failure of the BCS theory is a symptom of our 
inability to make accurate quantitative predictions in a 
strong coupling regime, rather than our having missed a 
qualitatively new type of order. Also, while this article 
will present a unified view of the important physics of 
the cuprate superconductors, it is not a comprehensive 
review, and it does not attempt to reflect the state of the 
field by representing the variety of viewpoints that have 
been taken elsewhere in the literature. 

The primary assertions of this article are as follows. 
At the lowest energy scales, the longest length scales, in 
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the absence of strong external perturbations, and at 'op- 
timal' carrier concentrations and above, all experimental 
indications are that the cuprate superconductors can in- 
deed be described in the framework of the BCS theory: 
the theory correctly captures the primary order parame- 
ter of the superconducting state, and the quantum num- 
bers of its elementary excitations. However, many ex- 
periments at lower doping concentrations and at shorter 
length scales require one or more additional order pa- 
rameters, either conventional (i.e. associated with the 
breaking of a sym metry o f the Hamiltonian) or 'topolog- 



cal' (see Section III.B.2j below). These order parame- 



ters are best understood and classified in terms of the 
physics of "Mott insulators," a topic which will be dis- 
cussed in greater detail below. The impo r tance of the 
Mott insulator was stressed by Anderson (1987). Our 



understanding of Mott insulators, and of their classifica- 
tion into categories with distinct physical properties has 
advanced greatly in the last decade, and a sharper ques- 
tion of experimental relevance is: which class of Mott 
insulators has its 'order' present in the cuprate super- 
conductors? As we shall discuss below, the evidence so 
far supports a class quite distinct from that implied in 



1991 



Anderson's proposal (Sachdev and Read 

How can the postulated additional order parame- 
ters be detected experimentally? In the simplest case, 
there could be long-range correlations in the new or- 
der in the ground state: this is apparently the case in 
La2-aSr,5Cu04 at low carrier concentrations, and we will 
describe recent experiments which have studied the in- 
terplay between the new order and superconductivity. 
However, the more common situation is that there are no 
long range correlations in any additional order parame- 
ter, but the 'fluctuating' order is nevertheless important 
in interpreting certain experiments. A powerful theoret- 
ical approach for obtaining semi-quantitative predictions 
in this regime of fluctuating order is provided by the the- 
ory of quantum phase transitions: imagine that we are 
free to tune parameters so that ultimately the new or- 
der does acquire long range correlations somewhere in 
a theoretical phase diagram. A quantum critical point 
will separate the phases with and without long-range or- 
der: identify this critical point and expand away from it 
towards the phase with fluctua ting order, which is the 
regime of experimenta l interest ( |Chubukov et at , 1994a 
jBachdev and~Yc] , gggg ); see Fig |l|. An illuminating dis- 



cussion of fluctuating order near quantum critical points 
(along with a thorough analysis of many recent experi- 
ments which has some overlap with our discussion here) 
has been provided recently by Kivelson et al. ( p002a ) . 



An especially important class of experiments involve 
perturbations which destroy the superconducting order of 
the BCS state locally (on the scale a few atomic spacings). 
Under such situations the theory outlined above predicts 
that the order of the Mott insulator is revealed in a halo 
surrounding the perturbation, and can, in principle, be 
directly detected in experiments. Perturbations of this 
type are Zn impurities substituting on the Cu sites, and 
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FIG. 1 Our theoretical strategy for describing the influence 
of a new order parameter in a BCS superconductor. Here g 
is some convenient coupling constant in the Hamiltonian, and 
we imagine that the superconductor of physical interest is a 
BCS superconductor with g > g c . Theoretically, it is useful to 
imagine that we can tune g to a value smaller than g c where 
there are long-range correlations in a new order parameter. 
Having identified and understood the quantum phase tran- 
sition at g = g c , we can expand away from it back towards 
the BCS superconductor (as indicated by the thick arrow) to 
understand the influence of quantum fluctuations of the new 
order parameter. This approach is most effective when the 
transition at g — g c is second order, and this will usually be 
assumed in our discussion. Note that the horizontal axis need 
not be the concentration of mobile carriers, and it may well be 
that the superconductor of physical interest does not exhibit 
the g < g c state at any carrier concentration. 



the vortices induced by an applied magnetic field. We 
shall discuss their physics below. 

To set the stage for confrontation between theory and 
experiment, we review some essential features of the BCS 
theory in Section ||, and introduce key concepts and or- 
der parameters in the theory of Mott insulators in Sec- 
We will combine these considerations in our 



in 



tion 

discussion of doped Mott insulators in Section |TV|, which 
will also include a survey of some experiments. A the- 
oretical phase diagram which encapsulates much of the 
physics discussed here appears in Section |y|, while Sec- 



VI concludes with a discussion on possible directions 



tion 
for future work 



II. BCS THEORY 



In BCS theory, superconductivity arises as an insta- 
bility of a metallic Fermi liquid. The latter state is an 
adiabatic continuation of the free electron model of a 
metal, in which all single particle states, labeled by the 
Bloch crystal momentum k, inside the /c-space Fermi sur- 
face are occupied by electrons, while those outside remain 
empty. With ct the creation operator for an electron 

with momentum k and spin projection a =||, a reason- 
able description of the Fermi liquid is provided by the 
free electron Hamiltonian 



#0 = H( £ fe ~ ^ C h C ka 



(1) 



ka 



where eg is the energy-momentum dispersion of the 
single-particle Bloch states and fi is the chemical poten- 
tial; the locus of points with eg = /j, defines the Fermi 
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surface. Changes in electron occupation numbers near 
the Fermi surface allow low energy processes which are 
responsible for the conduction properties of metals. 

BCS realized that an arbitrarily weak attractive inter- 
action between the electrons would induce the electrons 
near the Fermi surface to lower their e nergy b y bind - 
ing into pairs (known as Cooper pairs) (Cooper, 1956). 
BCS also proposed a mechanism for this attractive in- 
teraction: the exchange of a low energy phonon between 
two electrons, along with the rapid screening of the re- 
pulsive Coulomb interaction by the other electrons, leads 
to a residual attractive interaction near the Fermi sur- 
face. We regard this mechanism of electron pairing as a 
specific sidelight of BCS theory for good metals, and not 
an essential characterization of the BCS state. Indeed in 
liquid 3 He, the pairing is believed to arise from exchange 
of spin fluctuations ('paramagnons'), but the resulting 
superfluid state has many key similarities to the super- 
conducting metals. 

In the BCS ground state, the Cooper pairs undergo 
a process of condensation which is very closely related 
to the Bose-Einstein condensation of non-interacting 
bosons. Two well separated Cooper pairs obey bosonic 
statistics when adiabatically exchanged with each other, 
but their behavior is not simply that of point-like Bose 
particles when their internal wavefunctions overlap — the 
constituent electrons become important at these short 
scales; however it is the long distance bosonic charac- 
ter which is crucial to the appearance of a condensate of 
Cooper pairs. In the original Bose-Einstein theory, the 
zero momentum boson creation operator can be replaced 
by its c-number expectation value (due to the occupation 
of this state by a macroscopic number of bosons); simi- 
larly, the BCS state is characterized by the expectation 
value of the creation operator of a Cooper pair with zero 
center of mass momentum: 



lation: 



feT -H ki -fet 



oc A^ = Aq (cos k x — cos ky) . (2) 



The functional form of (g) in spin and fc-space carries 
information on the internal wavefunction of the two elec- 
trons forming a Cooper pair: we have displayed a spin- 
singlet pair with a d-wave orbital wavefunction on the 
square lattice, as is believed to be the case in the cuprates 
( |Scalapmc| , |1995b _ [Tsuei and Kirtley] , |2000| ). 

Along with (g) as the key characterization of the 
ground state, BCS theory also predicts the elementary 
excitations. These can be separated into two types: those 
associated with the motion of center of mass, R, of the 
Cooper pairs, and those in which a pair is broken. The 
center of mass motion (or superflow) of the Cooper pairs 
is associated with a slow variation in the phase of the 
pairing condensate A — > Aoe 1 ^': the superconducting 
ground state has 4>{R) — a constant independent of R 
(and thus long-range order in this phase variable), while 
a slow variation leads to an excitation with superflow. A 
vortex excitation is one in which this phase has a non- 
trivial winding, while the superflow has a non-zero circu- 



dR ■ V</> = 2Trn v 



(3) 



where n v is the integer-valued vorticity, and C is a con- 
tour enclosing the vortex core. A standard gauge invari- 
ance argument shows that each such vortex must carry 
a total magnetic flux of n v hc/(2e), where the 2e in the 
denominator represents the quantum of charge carried 
by the "bosons" in the condensate. Excitations which 
break Cooper pairs consist of multiple S = 1/2 fermionic 
quasiparticles with dispersion 



(4) 



and these reduce to the particle and hole excitations 
around the Fermi surface when Ao — > 0. 

All indications from experiments so far are that the 
cuprate superconductors do have a ground state char- 
acterized by (|^), and the elementary excitations listed 
above. However, BCS theory does make numerous other 
predictions which have been successfully and thoroughly 
tested in the low temperature superconductors. In par- 
ticular, an important prediction is that if an external 
perturbation succeeds in destroying superconductivity by 
sending |Ao| — > 0, then the parent Fermi surface, which 
was swallowed up by the Cooper instability, would reap- 
pear. This prediction is quite different from the perspec- 
tive discussed earlier, in which we argued for the appear- 
ance of a halo of order linked to the Mott insulator. 



III. MOTT INSULATORS 

The Bloch theory of metals also specified conditions 
under which crystalline materials can be insulating: if, 
after filling the lowest energy bands with electrons, all 
bands are either fully occupied or completely empty, then 
there is no Fermi surface, and the system is an insula- 
tor. However, some materials are insulators even though 
these conditions are not satisfied, and one-electron the- 
ory would predict partially filled bands: these are Mott 
insulators. Correlations in the motion of the electrons 
induced by their Coulomb interactions are crucial in pre- 
venting metallic conduction. 

One of the parent compounds of the cuprate super- 
conductors, La2Cu04, is a simple example of a Mott in- 
sulator. The lowest energy electronic excitations in this 
material reside on the Cu 3d x 2_ y 2 orbitals, which are lo- 
cated on the vertices of a square lattice. The crystal has 
a layered structure of stacked square lattices, with only a 
weak amplitude for electron hopping between successive 
layers. (We shall neglect the interlayer coupling and focus 
on the physics of a single square lattice in the remainder 
of this article.) After accounting for the ionization states 
of the other ions in La2Cu04, there turns out to be ex- 
actly one electron per unit cell available to occupy the 
Cu 3d 



X A —y z 



band. With two available spin states, this 
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FIG. 2 Motion of the two ferromagnetically aligned spins in 
(a) is prohibited by the Pauli principle. In contrast, the an- 
tiferromagnetically aligned spins at the top and bottom in 
(6) can access a high energy intermediate state (shown in the 
middle of (6)) and so undergo an exchange process. 



research in the last decade. We summarize the current 
understanding below. 

In keeping with the spirit of this article, we character- 
ize ground states of H by a number of distinct order pa- 
rameters. We only discuss states below which have long- 
range correlation in a single order parameter; in most 
cases, co-existence of multiple order parameters is also 
allowed ( patent s et~al\ , |1999| ; |Senthil and Fisher] , pOOOD , 
but we will ignore this complexity here. Our list of order 
parameters is not exhaustive, and we restrict our atten- 
tion to the most plausible candidates (in the author's 
opinion) for short-range Jij. 1 

Although our discussion below will refer mainly to 
Mott insulators, we will also mention ground states of 
non-insulating systems with mobile charge carriers: the 
order parameters we use to characterize Mott insulators 
can be applied more generally to other systems, and this 
will done in more detail in Section IV . 



band can accommodate two electrons per unit cell, and 
so is half-filled, and should have a metallic Fermi surface. 
Nevertheless, La2Cu04 is a very good insulator. The rea- 
son for this insulating behavior can be understood quite 
easily from a simple classical picture of electron motion 
in the presence of the Coulomb interactions. Classically, 
the ground state consists of one electron localized on each 
of the 3d x 2 _ y i orbitals: this state minimizes the repulsive 
Coulomb interaction energy. Any other state would have 
at least one orbital with two electrons, and one with no 
electrons: there is a large energetic penalty for placing 
two electrons so close to each other, and this prohibits 
motion of electrons across the lattice: hence the Mott 
insulator. 

Let us now look at the quantum theory of the Mott 
insulator more carefully. While charge fluctuations on 
each site are expensive, it appears that the spin of 
the electron can be rotated freely and independently 
on each site. However, in the quantum theory virtual 
charge fluctuations do occur, and these lead to residual 



"su per-ex change" interactions between the spins ( Andcr- 
sonT 195S). We represent the spin on the Cu site j by the 



5 = 1/2 spin operator Sj; the effective Hamiltonian that 
describes the spin dynamics then takes the form 



H - ^JijSi 



(5) 



where the are short-ranged exchange couplings and 
the ellipses represent possible multiple spin couplings, 
all of which preserve full SU(2) spin rotation invariance. 
Because the Pauli principle completely prohibits charge 
fluctuations between two sites if they have parallel spin 
electrons, while they are only suppressed by the Coulomb 
repulsion if they have opposite spins (see Fig ||), we ex- 
pect an antiferromagnetic sign > 0, so that nearby 
spins prefer opposite orientations. Classifying quantum 
ground states of models like (||) is a problem of consid- 
erable complexity, and has been the focus of extensive 



A. Magnetically ordered states 

Such states are obtained by examining H for the case 
of large spin S on each site: in this limit, the Sj can be 
taken as classical c-numbers, and these take a definite 
non-zero value in the ground state. More precisely, the 
SU(2) spin rotation symmetry of H is spontaneously bro- 
ken in the ground state by the non-zero values of (Sj), 
which are chosen to minimize the energy of H. We con- 
sider only states without a net ferromagnetic moment 
— 0), and this is expected because Jy > 0. The 
pattern of non-zero (Sj) can survive down to S = 1/2, 
and this is often found to be the case, although quan- 
tum fluctuations do significantly reduce the magnitude 

Of (Sj). 

An especially important class of magnetically ordered 
states 2 is characterized by a single ordering wavevector 



1 An order that has been much discussed in the literature, which 
we do riot, discuss here, is that, as sociated with the staggered flux 
state ( Affleck and Marston 198fr ) . and the rel ated algebraic spin 
liquid ( pantiicr and ^3^ 2001; |Wcrl , |2002a| ). The low energy 
theory nt these state s includes a rraple ^s I I M _ > gauge field, and 
it has been argued (Sachdev and Park, E003) that instantons, 
which are allowed because the underlying lattice scale theory 
has a compact gauge symmetry, always proliferate and render 
these u nstable toward s confining states (of the type discussed in 
Section III.B.l) in two spatial dimensions. However, states with 

in thr ee spatial dimensions 
2002b). 



a.pless TI(1) gauge field are 



(VIotrunich and Scnthil 



d are allowed i 
, |2002; ; [Wenj , f2> 



Magnetically ordered states with the values of (Sj) non-coplanar 
(i.e. three dimensional spin textures) are not included in this 
simple classification. Their physical properties are expected to 
be s imilar t o those of the non-collinear case discussed in Sec- 
tion [II. A. 2 in that quantum fluctuations of such a state lead to 
a paramagnet with topological order. However, this paramagnet 
is likely to have also a broken time-reversal symmetry. 
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FIG. 4 A state with non-collinear magnetic order on the 
square lattice defined by (^) and (Q) with wavevector K = 

(37T/4,7T). 



mobile carriers (Kivelson and Emery, 1996; Martin et al, 



20001; |Seibold et alj |1998|; [Tranquada et alj |1995t |Waki- 
': |\Y 



1999| ; |Whitc and Scalapinc| , |1998a| fl 
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FIG. 3 States with collinear magnetic order on a square lat- 
tice with unit lattice spacing and wavevectors (a) K = (n, 7r), 
(6) and (c) K — (37r/4, n). Shown are the values of (^|) on the 
square lattice sites Tj . A single unit cell is shown for the latter 
two states; they are crystographically inequivalent and have 
different reflection planes: in (b) the reflection planes are on 
certain sites, while in (c) they are at the midpoint between 
two sites. 



K: 

(Sj) = Ni cos (k ■ + N 2 sin (k ■ fj) (6) 

where fj is the spatial location of the site j, and N12 
are two fixed vectors in spin space. We list two key sub- 
categories of magnetically ordered Mott insulators which 
obey (p|): 



moto et al 



2001 



1999) — we can crudely view the mobile carriers as hav- 



ing induced an effective longer range exchange between 
the spins. Two examples of states with this value of K 
are shown in Fig |3[ a szie-centered state with N2 = in 
Fig^b, and a bond-centeied state with N2 = (\/2— l)Ni 
in Fig ||c. (The states have planes of reflection symme- 
try located on sites and the centers of bonds respectively, 
and so are crystallographically inequivalent. Also, these 
inequivalent classes are only present if the wavevector K 
is commensurate with the underlying lattice.) 



2. Non-collinear spins, Ni x N2 / 

Now the spin expectation values in (^|) lie in a plane 
in spin space, rather than along a single direction. For 
simplicity, we will only consider the simplest, and most 
common, case of non-collinearly ordered state, in which 



Ni-N 2 =0 ; N| = N 2 V0, 



(7) 



and then the values of (S , ) ma p out a circular spiral 
([ghraiman and Siggia , 1988 , 1989 ), as illustrated in Fig|J. 



1. Collinear spins, Ni x N 2 = 



In this situation, the mean values of the spins in (|6J) 
on all sites j are either parallel or anti-parallel to each 
other. The undoped insulator La2Cu04 is of this type 3 
with K = (tt, 7t); see Fig ||a. Insulating stat e s wit h 
static ho le s app e ared in Zaanen and Gunnarsson ( 1989 ), 
Machida| ( |1989| ), |5chufz| Ql989| ), and |Poilblanc and Rice 
(1989) with ordering wavevectors which move continu- 
ously away from (tt, tt). Another important illustrative 
example is the case K = (3tt/4, tt). Such a wavevector 
could be preferred in a Mott insulator by longer range Jy 
in (^|), but in practice it is found in a non- insulating state 
obtained by doping La2Cu04 with a suitable density of 



3 For this special value of K on the square lattice, and with the ori- 
gin of r co-ordinates on a lattice site, (a) is actually independent 
ofN 2 ) 



B. Paramagnetic states 

The other major class of states comprises those having 

(S,-) - 0, (8) 

and the ground state is a total spin singlet. 4 Loosely 
speaking each spin finds a partner, say Sj/, and the 
two pair up to form a singlet valence bond 



^(IT>m-'-U>ilT>i0 



(9) 



4 In a finite system with an even number of spins, the magnetically 
ordered ground state also has total spin zero. However, to obtain 
a state which breaks spin rotation symmetry as in (H) , it is neces- 
sary to mix in a large number of nearly degenerate states which 
carry non-zero total spin. The paramagnetic does not have such 
higher spin states available at low energy in a finite system. 
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FIG. 5 A crude variational wavefunction of a bond-ordered 
paramagnetic state. The true ground state will have fluctua- 
tions of the singlet bonds about the configuration shown here, 
but its pattern of lattice symmetry breaking will be retained. 
In other words, each bond represented by an ellipse above 
will have the same value of (Qai'Tj)}, and this value will be 
distinct from that associated with all other bonds. This pat- 
tern of symmetry breaking is represented more abstractly in 
Fig ft. ' 



Of course, there are many other choices for the part- 
ner of spin Sj, and in the Feynman path integral pic- 
ture we imagine that the pairing configuration fluctuates 
in quantum imaginary time; t h is is th e 'resonating va- 



lence bond' picture of Pauling (1949), Fazekas and An- 
derson] (1974), and Anderson (1987). However, there is 
a great of structure and information contained in the 
manner in which this fluctuatio n takes place, and re- 
search (Chubukov et a/. , 1994b , 3; Read and Sachdev 



1991 ; Sachdev and Rcaq 7 1991 ) dilineating this structure 



has led to the following classification of paramagnetic 
Mott insulators. 



1. Bond-ordered states: confined spinons 

This class of states can be easily understood by the 
caricature of its wavefunction shown in Fig J5f here each 
spin has chosen its valence bond partner in a regular 
manner, so that there is a long-range 'crystalline' order 
in the arrangement of valence bonds. This ordering of 
bonds clearly breaks the square lattice space group sym- 
metries under which the Hamiltonian is invariant. Of 
course, the actual wavefunction for any realistic Hamil- 
tonian will have fluctuations in its valence bond configu- 
ration, but the pattern of lattice symmetry breaking im- 
plied by Fig will be retained in the true bond-ordered 
ground state. We can make this precise by examining 
observables which are insensitive to the electron spin di- 
rection: the simplest such observables we can construct 
from the low energy degrees of freedom of the Mott in- 
sulator are bond variables, which are a measure of the 
exchange energy between two spins: 



(a) 




nu:: 



- - - 



v ■ - - 



v 



(C) 



□ i:::: 



-> - 



(d) 



FIG. 6 Pattern of the bond variables {Qa(rj)}, for f a a 
nearest-neighbor vector, in a number of paramagnetic states 
with (Sj) = 0. For each state, the values of {Qa(fj)) are equal 
on bonds represented by the same type of line, and unequal 
otherwise. The number of distinct values of (Q a {fj)) are (a ) 
3, (b ) 2, (c ) 5, and (d ) 5. The unit cells of the ground states 
have sizes (a ) 2 x 1, (b ) 2 x 2, (c ) 4 x 1, and (d) 4x4. 



a measure of the charge density 5 on site rj , and so this 
special case of ( |l0| ) measures the "charge order." 

The state introduced in Fig [| can be characterized by 
the pattern of values of (<3 a (r}j) with r a a nearest neigh- 



)a(fj) = Sj ■ S 



j+a- 



(10) 



Here a denotes displacement by the spatial vector r a , 
and so the spins above are at the spatial locations fj and 
fj + f a . We will mainly consider bond order with r a =/= 0, 
but note that the on-site variable Qo(fj), with r a = 0, is 



5 By (|10|), Qo(fj) = S?. A site with a spin has S 2 - = 3/4, while a 
site with a hole has = 0, and we assume that doubly occupied 

sites are very rare. Thus S?, and hence Qo(fj) is seen to be 
linearly related to the charge density on site j. 
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bor vector, as shown in Figga: notice there are 3 distinct 
values of (Q a (fj)) and symmetries of the states in Figs || 
and ||a are identical. While these 3 values are quite dif- 
ferent in the trial state in Fig their values in the actual 
ground state may be quite close to each other: it is only 
required that they not be exactly equal. 



„ — „ , 

appeared in some theories (Altman and Auerbach 


2002 


Dombrc and Kotliai 


M: 


Read and Sachdcv 




1989a 


Sachdev and Read, 199(i|), 


is shown in Fig 1(3 


d: here 



the bonds have a plaquette-like arrangement rather than 
columnar, but, as we shall discuss below, the physical 
properties of all the states in Fig ^ are quite similar to 
each other: 



We can also consider patterns ot bond order with larger 
unit cells, and two important structures which have ap- 
peared in theories o f doped Mott insulators ( Vojta , 2002 ; 
Vbjta and Sachdev, 1999) are shown in Figs ||c and 
(related bond orders also appear in studies of quasi-one- 
dime nsional models appropriate to organ i c sup erconduc- 
tors (|Clay et al{ |2002| ; pVIazumdar et al\ , |200C| )). Again, 
as in Section III.A.l, such states could, in principle, also 
appear in Mott insulators with longer-range exchange in 
(||). An interesting property of these states is that, un- 
like the states in Fig ||a and ||b, not all sites are crystal- 
lographically equivalent. This means that on-site spin- 
singlet observables, such as the site charge density, will 
also have a spatial modulation from site to site. A sub- 
tlet y is that the Hamiltonian acts on a Hilbort space 



of 51 — 1/2 spins on every site, and so the charge den - 



sity on each site is fixed at unity. However, it must 
be remembered that (^|) is an effective model derived 
from an underlying Hamiltonian which Hoes allow vir- 



tual ! charge fluctuations, and the site charge modulations 
in the states of Figs ^c and ||d will appear when it is 
properly computed in terms of the microscopic degrees 
of freedom. At the same time, this argument also makes 
it clear that any such modulation is suppressed by the 
repulsive Coulomb energy, and could well be difficult to 
observe, even in the doped antiferromagnet. So the on- 
site variable, (Qo(fj)) = (S|), will have a weak modula- 
tion in the states of Fig ||c and |6|d when computed in the 
full Hilbert space of the model with charge fluctuations. 
Note, however, that the modulation in bond orders asso- 
ciated with Q a {fj), with f a ^ 0, need not be small in the 
states in Fig ^, as such modulations are not suppressed 
as effectively by the Coulomb interactions. 

The physical mechanism inducing bond-ordered states 
such as those in Fig || is illustrated in the cartoon pic- 
tures in Fig |?]. More detailed computations rely on a 
semiclassical theory of q uantum fluctuations near a mag- 
netically ordered state (Read and Sachdcv, 199C). Re- 
markably, very closely related theories also appear from 
a very different starting po i nt — from duality mappin gs 
( Fradkin and Kivclson , 1990|; Read and Sachdev , 1990|) o f 
"quantum dimcr models" (Rokhsar and Kivelson, 1988) 
of the paramagnetic state. These computations show 
that spontaneous bond order invariably appears in the 




FIG. 7 Bond order induced by quantum fluctua tions. Va- 



lence bonds gain energy by "resonating" in pairs (Anderson 



1987; Fazekas and Anderson 



197-1 



Pauling 



1949 



Rokhsar 



and Kivelson, 1988| ); shown are resonances around the pla- 
quette (i.e. square loop) marked with a star. For the regular 
bond-ordered configuration of valence bonds in (a), such res- 
onance can occur not only around the plaquette marked with 
a star, but around five additional plaquettes. In contrast, 
in (b), such a resonance is possible only around the plaque- 
tte marked with a star. This additional quantum "entropy" 
associated with (a) selects regular bond order in the ground 
state. More sophisticated considerations (which also allow va- 
lence bonds that do not connect nearest neighbor sites) show 
t hat this mechanism is es p ecially effective in two d imensions 
(|Read and Sachdevj, [1990] fcachdev and Park[ §00§). 



ground state in systems with collinear spin correlations in 
two spatia l dime nsions (Read and Sachdev, 1990; Sachdev 
and Park , [2002 ). We will have more to say about this 



connection between and bond and collinear spin order in 



Section 1II.C.1 



We also mention here the "nematic" states of Kivcl 



son et al.\ (|1998[ ) in the doped Mott insulator. These 



can also be characterized by the bond order variables in 
(|l0|). The symmetry of translations with respect to fj 
is not broken in such states, but the values of (Q a {fj)) 
for symmetry-related values of f a become unequal. For 
example, (Q a (fj)) has distinct values for f a = (1,0) and 
(0, 1). Such states also appear in certain insulating anti- 
ferromagnets ( |Read and Sachdev] , |1989a| Jb|, |1990| ). 

It also interesting to note here that the bond order vari- 
ables Q a (fj ) also have spatial modulations in some of the 
magnetically o rdere d states considered in Section III. A 
( p^achar et al , 1998). It is clear from (|h]) that any bro- 
ken lattice symmetry in the spin-rotation invariant quan- 
tity (Sj) ■ (Sj +a ) will generate a corresponding broken 
symmetry in the bond variable {Q a {fj))- Evaluating the 
former using (||) we can deduce the following: (i) the 
K = (tt, 7r) state in Fig ||a and the spiral state in Fig |] 
have {Qaifj)) independent of fj, and hence no bond or- 
der; (ii) the bond-centered magnetically ordered state in 
Fig ||c has precisely the same pattern of bond order as 
the paramagnetic state in Fig ||c; (m) the site-centered 
magnetically ordered state in Fig ||b has bond order with 
(Qa^)) fj .-dependent, but with a pattern distinct from 
any shown here — this pattern of bond order is in principle 
also allowed for paramagnetic states, but has so far not 
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FIG. 8 Linear confining potential between two neutral S = 
1/2 spinons in a bond-ordered state. The line of valence bonds 
with dashed lines is out of alignment with the global bond 
order, and it costs a finite energy per unit length. 



been found to be stable in various studies. Finally, note 
that in (ii) and (Hi) the period of the bond order (four) 
is half that of the spin modulation (eight)-this is easily 
seen to be a general relationship following from the corre- 
spondence (Q a (fj)) ~ (Sj) • (Sj-+ ) + •■■ in magnetically 
ordered states, which with (^3j) implies an j^— dependent 
modulation of the bond order with wavevector 2K. It is 
worth reiterating here that this last relationship should 
not be taken to imply that there are no modulations in 
(Qai'Tj)) when (S„) = 0: there can indeed be bond mod- 
ulations in a paramagnet, as discussed in the other para- 
graphs of this subsection, and as is already clear from 
the simple wavefunction in Fig || — these will be impor- 
tant later for physical applications. 

We continue our exposition of paramagnetic bond- 
ordered states by describing excitations with non-zero 
spin. These can be understood simply by the analog of 
cartoon wavefunction pictures drawn in Fig |5| To create 
free spins we have to break at least one valence bond, 
and this initially creates two unpaired, neutral, 5=1/2 
degrees of freedom (the "spinons"). We can ask if the 
spinons can be moved away from each other out to in- 
finity, thus creating two neutral 5 = 1/2 quasiparticlc 
excitations. As illustrated in Fig || this is not the case: 
connecting the two spinons is a line of defect valence 
bontls which is not properly aligned with the global bond 
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FIG. 9 Cartoon wavefunction for 2 static spinless Zn impuri- 
ties in a confining, bond-ordered state. It we attempt to con- 
struct a wavefunction only using singlet valence bonds, then 
just as in Fig ^, there will be defect line of singlet bonds which 
are not aligned with the global bond order, which will cost a 
finite energy per unit length. When the two Zn impurities 
are sufficiently far apart, it will pay to restore the bond order 
in between the impurities, at the price of unpaired S = 1/2 
moments, one near each impurity. 



of removing a 5 = 1/2 spin from a fixed site j in (j^). 
Experimentally, this can be conveniently done by substi- 
tuting a spinless Zn ++ ion in place of an 5 = 1/2 Cu ion. 
The main physical effect can be understood from the car- 
toon wavefunction in Fig |^: it is convenient to imagine 
placing 2 Zn impurities, and then moving them apart out 
to infinity to deduce the physics in the vicinity of a sin- 
gle impurity. As in our discussion above for spinons, note 
that there will initially be a line of defect valence bonds 
connecting the two Zn impurities, but it will eventually 
pay to annihilate this defect line by creating two spinons 
and binding each to a Zn impurity. Thus each Zn im- 
purity confines a free 5=1/2 spinon in its vicinity, and 



this can be detected in experiments (Finkelstcin et al 
199fj|). 6 



2. Topological order: free spinons 

This type of paramagnet is the "resonating valence 
bond" ( RVB) state (Anderson, 1987; Baskaran and An 



dersorj, 1988: Fazckas and Anderson , 1974; Kivclsou 



orde jr, and these detects have a finite energy cost per unit e t g |1987| ; |Mocssncr and Sondhi| , |200T| gaujjnj |1949| ) 



length" This linearly increasing potential is quite anal- 
ogous to that between a quark and an anti-quark in a 
meson, and the spinons (quarks) are therefore perma- 
nently confined (Read and Sachdcv, 1989b). Moving two 
spinons apart from each other will eventually force the 
breaking of the defect line by the creation of another 
pair of spinons. The only stable excitation with nonzero 
spin therefore consists of a pair of spinons and carries 
spin 5=1. We will refer to this quasiparticle as a spin 
exciton as its quantum numbers and observable charac- 
teristics are quite similar to spin excitons found in semi- 
conductors and metals. The spin exciton is clearly the 
analog of a meson consisting of a quark and anti-quark 
pair. 

A similar reasoning can be used to understand the in- 
fluence of static spinless impurities i. e. the consequences 



in which the singlet pairings fluctuate in a liquid-like con- 
figuration, in contrast to the crystalline arrangement in 



In principle the Zn impurity could also bind an electron (with 
or without a spinon) but this is suppressed by the charge gap 
in a Mott insulator. Later, in Section IV. C when we consider 



Zn impurities in d-wave superconductors, a related phenomenon 
appears in the fo rm nf the Kondri effect. 
7 In recent years, Anderson (2002) has extended the RVB con- 



cept to apply to doped Mott insulators at temperatures above 
T c . This extension is not in consonance with the classification 
of the present article. The topological order discussed in this 
subsection can only be defined at T = in two spatial dimen- 
sions. The description at T > T c requires solution of a problem 
of quantitative difficultly, and with incoherent excitations, but 
without sharp distinctions between different states. 



9 




FIG. 10 Topological order in a resonating valence bond state. 
Shown is one component of the wavefunction, with a partic- 
ular pairing of the spins into local singlets: the actual wave- 
function is a superposition over a very large number of such 
pairing configurations. The number of valence bonds cutting 
the dashed line is an invariant modulo 2 over these pairing 
configurations, as shown by tne following simple argument. 



of this Z2 gauge theory to Section III.C.2| (see especially 
Fig |ll|): here, we show that such a gauge theory has 
similar topological properties. In a system with periodic 
boundary conditions (with the topology of a torus) , the 
Zi gauge theory has different sectors depending upon 
whether there is a Z2 flux piercing a ny of the holes of the 
torus (following [Benthil and Fisher (2000), this Z2 flux is 
now commonly referred to as a "vison"). In the valence 
bond picture discussed in the previous paragraph, a vi- 
son changes the sign associated with every valence bond 
cutting a line traversing the system in the vison direc- 
tion (the dashed line in Fig [)]); in other words, the even 
and odd valence bond sectors mentioned above now have 
their relative signs in the wavefunction changed. 

In addition to appearing in the holes of the torus, the 
vison can also appear as a singlet excitation with i n the 



bulk (Kivclson, 198S 



Any leaiiaiigemeiil of tire valence bonds can be leached by 

repeated application of an elementary rearrangement between 
4 spins: (1, 2)(3, 4) — » (1,3)(2,4) (here (i,j) denotes a singlet 
bond between Si and Sj). So it is sufficient to check this con- 
servation law for 4 spins: this is done easily by explicitly con- 
sidering all different possibilities among spins 1,2,3,4 residing 
to the left /right of the dashed line. If the system has periodic 
boundary conditions along the horizontal direction, then this 
conservation law is violated, but only by rearrangements as- 
sociated with loops which circumnavigate the systems; these 
only occur with a probability which becomes exponentially 
small as the circumference of the system increases. 



Fig |5|. Despite the apparent 'disorder' in the valence 
bond configuration in the ground state, there is actually 



this type of Mott insulator ( 


Boncstccl|, 1986 ; Kivclson 


1989 


; |Read and Chakraborty 


, 1989; Read and Sachdcv 


1991 


; Rokhsar and KivelsoE 


> 


1988]; ThoulesE, 


1987; 


Wen, 


1991 


) , and which plays an important role in determining 



its excitation spectrum. The reader can see this in the 
context of the cartoon picture shown in Fig Count 
the number of singlet valence bonds cutting the dashed 
line in this figure: this number will clearly depend upon 
the particular valence bond configuration chosen from 
the many present in the ground state, and one such is 
shown in Fig [Io|. However, as argued in the figure cap- 
tion, the number of bonds cutting the dashed line is con- 
served modulo 2 between any two configurations which 
differ only local rearrangements of valence bonds: the 
quantum number associated with this conservation is the 
topological order in the ground state. 

A convenient and powerful description of this topolog- 
ical order is provided by an effective model of the singlet 



sector formulated as Z2 gauge t heory ( Read and Sachdcv 



1991|; Bachdev and Rcadj |1991| |Scnthil and Fisher| , |2000 
1991 jy We postpone a self-contained derivation 



Read and Chakraborty , 1989; Read 



and Sachdev| , |1991j Benthil and Fisher , 2000 ) . It is now 
a vortex excitation in the Z2 gauge theory, that requires 
a finite energy for its creation. We will see below in 



Section [II. C. 2 that there is an alternative, and physi- 
cally revealing, interpretation of this vortex excitation in 
terms of the order parameters used earlier to characterize 
the magnetically ordered state, and that the topological 
order is intimately connected to the vison energy gap. 

Finally, we can describe the spin-carrying excitations 
of this topologically ordered state using the crude, but 
instructive, methods used in Section [II.B.l. As there 



is no particular bond order associated with the ground 
state, the spinons have no confining force between them, 
and are perfectly free to travel throughout the system as 
independent neutral S — 1/2 quasiparticles. Similarly, 
there is no confining force between Zn impurities and the 
spinons, and so it is not required tha t an 5* = 1/2 m oment 



be present near each Zn impurity (Fendley et al 
Bachdev and Vojta| , |2000| ). 



2002 



C. Connections between magnetically ordered and 
paramagnetic states 

A central ingredient in the reasoning of this article 
is the claim that there is an intimate connection be- 



tween the magnetically ordered states in Section III. A 
and a corresponding paramagnetic state in Se ction [TLB . 



In particular, the collinear states of Section III.A.l 



linked to the bond-ordered states in Section III.B.l, while 



Readers not familiar with Z2 gauge theories may understand 



them by analogy to electromagnetism. The latter is a U(l) gauge 
theory in which the physics is invariant under the transforma- 
tion 2 — > e 1( ^z, A^ — > Ajj — d^<f> where z is some matter field, 
A/j, is a gauge field, and is an arbitrary spacetime-dependent 
field which generates the gauge transformation. Similarly, in a 
Z2 gauge theory, matter fields transform as z — > rjz, where r\ is 
a spacetime-dependent field which generates the gauge transfor- 
mation, but is now allowed to take only the values r\ = ±1. The 
Z2 gauge field 0Y7 resides on the links of a lattice, and transforms 
as Uij rjiCTij-qj. 
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the non-collinear states of Section III. A. 2 are link 



3d to 
The 
llow- 
ithor 



the topologically ordered states of Section [II. B. 2 
reader will find a more technical discussion of the fc 
i ng issues in a companion review article by the a' 
(|Sachde4 [20031 ) . 

Before describing these links in the following sdbsec- 
tions, we discuss the meaning of the "connectedness 
two states. The magnetically ordered phases are 
acterized by simpl e orde r parameters that we hav 
cussed in Section [II. A. Now imagine a second- 



quantum phase transition in which the magnetic 
range order is lost, and we reach a state with 
ating magnetic correlations, which is ultimately a 
tionally invariant, spin-singlet paramagnet at the lo 
length scales. We will review arguments below 



of 
char- 
dis- 
arder 
long- 
u- 
rota- 
agest 
■\ ifhich 



fl ncti 



show that t his "quantum disordered" state ( Chakra|varty 
et aj], |1989| ) is characterized by the order parameter i >f the 
connected paramagnetic state i.e. fluctuating col. incar 
magnetic order leads to hour] order, while fluctuating 



non-lcollincar magnetic order can lead to topologic 



iL 



der. So two connected states are generically proximite to 
each other, without an intervening first order transition 
in a generalized phase diagram drawn as a functisn of 
the couplings present in the Hamiltonian. 



1. Collinear spins and bond order 



It should be clear from Section [II. A. 1 that coljincar 

The 
l by 



N 



tances 



fluct 



spin states are characterized by a single vector Ni 
second vector N2 is pinned to a value parallel to 
some short distance physics, and at long dist 
may consider a theory of the fluctuations of Ni ; 
In a phase with magnetic order, the dominant spin- 
fluctuations occur in configurations with a fixed nor 
value of |Ni|. In the transition to a non-magnetic 
the mean value of |Ni| will decrease, until the 
tions of Ni occur about Ni — in a paramagnetic 
There are 3 normal modes in this fluctuation spec 
corresponding to the 3 directions in spin space, an 
resultant is an S — 1 gapped quasiparticle excitatipn 
the paramagnetic state. This we can easily identify 
S = 1 spin exciton of the bond-ordered state: this 
tification is evidence supporting our claimed connc 
between the states of Section III.A.l and 



[II.B.l 



Further evidence is provided by detailed computations 
which show the appearance of bond order in the regime 
where Ni fluctuations have lost their long-range ( rder. 
We have already seen a simple example of this abc ve in 
that the magnetically ordered state in Fig ||c ahead r had 
the bond order of the paramagnetic state in Fig ^c it is 
completely natural for the bond order in the magnet cally 
ordered phase in Fig |^c to persist across a transiti an in 
which spin rotation invariance is restored, and this con- 
nects it to the state in Fig ^c. A non-trivial exa mple 
of a related connection is that between the K — tt, 7t) 
Neel state in Fig ||a, and the paramagnetic bond-or lered 
states in Figs @a and ||b, which was established by [lead 



and Sachdev 



£21= 



we 
lone, 
wave 
zero 
p|hase, 
ua- 
p hase. 
rum, 
i the 
in 
the 
iden- 
ction 



i & 1 



Sachdev and 



Park ( 2002 ): Berry phases associated with 
of the lattice spins were shown, after a 
mapflimg, to induce bond order in the phase in 



the precessiojn 
duality 

which long-r;jige order in Ni was lost. 



2. Non-colline;ir spins and topological order 



The first ai gument of Section III.C.l , when generalized 



to non-collin( ar 
subtle characterization 
phase. 

Recall fron l 
nctic phase is 
length, vectors 
two vectors 
ber of real 
state to 4. Tliere 



Section [II.A.2| that the non-collinear mag- 
characterized by two orthogonal, and equal 
N 12 . It takes 6 real numbers to specify 
>ut the 2 constraints in (Q) reduce the num- 
rameters required to specify the ordered 
is a useful parameterization ( Chubukov 
et al\ 1994b , ]) which explicitly solves the constraints (Q) 
Ni : 2 in terms of 2 complex numbers Zj-, z\ 
c uivalent to the required 4 real numbers): 



I ar; 



by expressin 
(which are e 



' wii h 



( 



It can also b 
like an S = 
of dealing 
we can expre 

cial price we 
that the 
the spinors 
the same non 
change the 
spacetime wi 
tive action fo 
invariance 
logically or> 
where we ha< . 
theory. 

In the ma; 
pect dominaJit 



non-zero value 



stramt |Zf 
sphere in a f 
ordered groupd 
components 
posite points 
and {—z-\, 
order 

has crucial 
tations that 
dered and 
order parameter 
(|19790 ) allows 
homotopy gr< >up 



I) 

parameter 



tie 



(|l989b|), [lead and Sachdevj ( |1990D , and 



spins, leads quite simply to a surprisingly 
of the associated paramagnetic 



Ni + iN 2 



7 2 - 7 2 

■ 2 2 

2^1 



(11) 



checked from (pd| ) that (zf, z\) transforms 
/2 spinor under spin rotations. So instead 
a constrained theory of Ni,2 fluctuations, 
>s the theory in terms of the complex spinor 
whiqh is free of constraints. There is one cru- 
have to pay for this simplification: notice 
par^metrization ( |ll| ) is double-valued and that 
Zf,zjJ and (— Zf,— z{) both correspond to 
collinearly ordered state. Indeed, we can 
of z independently at different points in 
hout changing the physics, and so any effec- 
the (zf, z\) spinor must be obey a Zi gauge 
Here is our first connection with t he topo- 
d$red paramagnetic state of Section III.B.2 , 
also discussed a description by a Z 2 gauge 



sign 



netically ordered non-collinear state we ex- 
rotational fluctuations about some fixed 
Nf = Nl = (|z T | 2 + |z x | 2 ) 2 . The con- 
zj 2 = constant defines the surface of a 
)ur-dimensional space (S3) of magnetically 
states defined by the real and imaginary 
)f z-f, z\. However, we need to identify op- 
on the sphere with each other, as (z-\,z{) 
) are equivalent states: this identifies the 
space with S3/Z2. This quotient form 
consequences for the topological defect exci- 
ire permitted in both the magnetically or- 
paramagnetic phases. In particular, the 
space (see the review article by Mcrmin 



stable Z 2 vortices associated with the first 
m (S3/Z2) — Z2: upon encircling such a 
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FIG. 11 A vison flSenthil and Fisher| , (200C| ). On the left we 
show a circular path in real space; this path could be entirely 
within the bulk of the system (in which case it defines a local 
vison excitation) or it encircles the entire system, which obeys 
periodic boundary conditions (so that now it defines a global 
topological excitation). On the right is the space of magneti- 
cally ordered states represented by the complex spinor (zj , Zj) 
up to an overall sign. As we traverse the real space circle, the 
path in order parameter space connects polar opposite points 
on S3 (A and B), which are physically indistinguishable. A 
key point is that this vison excitation can be defined even in 
a state in which magnetic order is lost: the path on the right 
will fluctuate all over the sphere in quantum imaginary time, 
as will the location of the points A and B, but A and B will 
remain polar opposites. 



vortex, we traverse a path in the order parameter space 
from (z-|-, Zi) to (—z^,—z{), as shown in Fig[ll]. As argued 
in the caption, a fundamental point is that such vortices 
can be defined as sensible excitations even in the para- 
magnetic phase, where (zj, z{) is strongly fluctuating in 
quantum imaginary time: upon encircling the vortex, the 
path in order parameter space will also strongly fluctu- 
ate, but will always connect polar opposite points on S3. 
We identify these paramagnetic vortices with the visons 
of Section III.B.2, thus firmly establishing a connection 
between non-collinear magnetic order and the topologi- 
cally ordered paramagnet. 



Finally, we wish to consider a Z2 gau ge theory in which 



magnetic order is lost continuo usly (Chubukov et al. 



1994c| ; |Read and Sachdev| , |1991| ), and we obtain a para- 
magnetic phase in which the spinor (zf,Z\) fluctuates 
about 0. A ped agogical descript i on of such a theory 
was p rovid ed by Lammert et al\ (1993) and Lammert 
et a{] ( |199E ) in an entirely different context: they consid- 
ered thermal phase transitions in a nematic liquid crystal, 
with order parameter S2/Z2, in three spatial dimensions. 
However their results can be transposed to the quantum 
phase transition in two spatial and one imaginary time 
dimension of interest here, with the primary change being 
in the order parameter space from S2/Z2 to S3/Z2: this 
change is only expected to modify uninteresting numeri- 
cal factors in the phase diagram, as the global topologies 
of the two space s are the same. As sh own by Lammert 
et a{] ( |1993 ) and Lammert et al. ( 1995| ), the magnetically 
ordered state (with states labeled by points in £2,3/^2) 
does indeed undergo a continuous phase transition to a 
paramagnetic state in which spin rotation invariance is 
restored and a topological order is present. This topo- 



logical order arises because the Z2 visons discussed in 
Fig [ll] do not proliferate in the paramagnetic state; in 
this sense, the topological order here is similar to the 
topological order in the low temperature phase of the 
classical XY model in 2 dimensions, where point vortices 
are suppressed below the Kosterlitz-Thouless transition 



(FhoulesE, 1998). We can also connect th e nonpr olifcra- 



tion of visons to our discussion in Section [H.B.2| > where 



we noted that there was an excitation gap tow ards the 
creation of Z2 visons flScnthil and Fishcij |200C| ). Indeed, 
an explicit connection between the topological order be- 
ing discussed here and the topological order noted in the 
captio n to Fig | Io| was esta b lished by Read and Sachdcv 



(1991), Sachdcv and Read (1991), and Chubukov et al. 



(|1994rJ). 

Moreover, without the proliferation of visons in the 
ground state, the (z-\,zi) configurations can be defined 
as single-valued configurations throughout the sample. 
Normal- mode oscillations of (z-\,z±) about zero can now 
be identified as a neutral S = 1/2 particle. This is clearly 
related to the spinon excitation of Section III.B.2 ; this is 
our final confirmation of the intimate connection between 



the non-collinear magnetic states of Section III. A. 2 and 
the topologically ordered states of Section III.B.2. 



This is a good point to mention, in passing, recent 
neutron scatte r ing ev idence for a RVB state in CS2CUCI4 
( |Coldca et al , 2001); the measurements also show non- 
collinear spin correlations, consistent with the connec- 
tions being drawn here. 



IV. ORDER IN STATES PROXIMATE TO MOTT 
INSULATORS 

We are now ready to discuss the central issue of order 
parameters characterizing the cuprate superconductors. 
These superconductors are obtained by introducing mo- 
bile charge carriers into the Mott insulator of the square 
lattice of Cu ions that was discussed at the beginning 
of Section III. The charge carriers are introduced by 
substitutional doping. For instance, in the compound 



La2-aSr ( 5Cu04, each trivalent La 3+ ion replaced by a di- 
valent Sr 2+ ion causes one hole to appear in the Mott 
insulator of Cu ions: the concentration of these holes is 
6 per square lattice site. 

For large enough S, theory and experiment both in- 
dicate that such a doped Mott insulator is a d-wave 
superconductor characterized by the pairing amplitude 
(^). The reader can gain an intuitive (but quite crude 
and incomplete) understanding of this by the similar- 
ity between the real-space, short-range pair in (^J) and 
the momentum-space, long-range pairing in (|^). The 
undoped Mott insulator already has electrons paired 
into singlet valence bonds, as in (|9|), but the repulsive 
Coulomb energy of the Mott insulator prevents motion 
of the charge associated with this pair of ele ctrons . It 
should be clear from our discussion in Section III.B that 



this singlet pairing is complete in the paramagnetic Mott 
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insulators, but we can also expect a partial pairing in the 
magnetically ordered states. Upon introducing holes into 
the Mott insulator, it becomes possible to move charges 
around without any additional Coulomb energy cost, and 
so the static valence bond pairs in (^) transmute into the 
mobile Cooper pairs in (||); the condensation of these 
pairs leads to superconductivity. Note that this discus- 
sion is concerned with the nature of the ground state 
wavefunction, and we are not implying a "mechanism" 
for the formation of Cooper pairs. 

The discussion in the previous sections has laid the 
groundwork for a more precise characterization of this 
superconductor using the correlations of various order 
parameters, and of their interplay with each other. The 
proximity of the Mott insulator indicates that the Cooper 
pairs should be considered descendants of the real-space, 
short range pairs in (||), and this clearly demands that 
all the magnetic, bond and topological order parameters 
discussed in Section remain viable candidates for the 
doped Mott insulator. The motion of charge carriers al- 
lows for additional order parameters, and the most im- 
portant of these is clearly the superconducting order of 
the BCS state noted below (||) in Section El In principle, 
it is also possible to obtain new order parameters which 
are characteristic of neither the BCS state nor a Mott 
insulator, but we such order parameters shall not be dis- 
cu ssed her e (discussions of one such order may be found 
i n [Hsu et_ 



( 2003 ), Chakravarty et al. ( 2001 ), and Bchollwock et al 
(2003)) 



(|1991|), |Wcn and Leej ( |1996D, |Lee and Sha 



The arsenal of order parameters associated with Mott 
insulators and the BCS state permits a very wide va- 
riety of possible phases of doped Mott insulators, and 
of quantum phase transitions between them. Further 
progress requires experimental guidance, but we claim 
that valuable input is also obtained from the theoretical 



connections sketched in Section III.C 



1991 



The simplest line of reasoning ( Bachdev and Reac] 



uses the fact that the undoped Mott insulator 
La2Cu04 has collinear magnetic order as sketch ed in 
Fig ||a. The arguments above and those in Section III.C 
then imply that the doped Mott insulator should be 
characterized by the collinear magnetic order of Sec- 



tion 1II.A.1, the bond order of Section III.B.l, along with 
the phase order of BCS theory. This still permits a large 
variety o f phase diagrams, and some of these were ex- 
plored in Bachdev and Read ( 1991 ), Vojta and Sachdcv 
( |1999|) , |Vqjta et al\ (|2000aD , and |Vojta| ( |2002|) , with ck7 
tailed results on the evolution of bond order and super- 
conductivity with increasing doping. However, this rea- 
soning excludes phases as sociated with the non-collincar 
magnetic orde r of Sec tion III. A. 2 and the topological or- 
der of Section pILR2 . 



Some support for this line o f reasoning came f r om th e 



breakthrough experiments of Tranquada et al. _[1995) 



Tranquada et al. ( 1996| ), and Tranquada et al. (1997) on 
L^-y-aNdySraC^ for hole concentrations near 5 = 1/8: 
they observed static, collinear magnetic order near the 



wavevectors K = (3ir/4, it) shown in Figs ||b,c, which co- 
existed microscopically 9 with superconductivity for most 
5. They also observed modulations in the bond order 
Q a (rj) (Eqn (JToj) ) at the expected wavevector, 2K. The 
experimentalists interpreted their observations in terms 
of modulations of the site charge density — proportional 
to Qo(rj) — but the existing data actually do not discrim- 
inate between the differ ent poss ible values of r a . As we 
noted earlier in Section III.B.l, the physical considera- 



tions of the present article suggest that the modulation 
may be stronger with r a ^ 0. (The existing data also can- 
not distinguish between the magnetic orders in Fig [p3 



(site-centered) and Fig 
the bond orders in Fig 



(bond-centered), or between 
(orthorhombic symmetry) and 
Fig |6|d (tetragonal symmetry).) We also mention here the 
differen t physical considerations in the early theoretica l 
work o f [Zaancn and Gunnarsson ( 1989 ), Machida (1989). 
5chulz| ( jl989| ), and |Poilblanc and Ricej ( |1989| ) which led 
to insulating states with collinear magnetic order with 
wavevector K ^ (tt 7 it) driven by a large site-charge den- 
sity modulation in the domain walls of holes. 

The following subsections discuss a number of recent 
experiments which explore the interplay between the or- 
der parameters we have introduced here. We argue that 
all of these experiments support the proposal that the 
cuprate superconductors are characterized by i nterplay 



between the collinear magnetic order of Section III.A.l 



the bond order of Sec tion III.B.l (these are connected 
as discussed in Section III.C.l), and the superconducting 
order of BCS theory. 



A. Tuning order by means of a magnetic field 

In Section |, we identified a valuable theoretical tool 
for the study of systems with multiple order parameters: 
use a coupling g to tune the relative weights of static or 
fluctuating order parameter correlations in the ground 
state. Is such a coupling available experimentally ? One 
choice is the hole concentration, 5, and we can assume 
here that g increases monotonically with S. However, 5 
is often difficult to vary continuously, and it may be that 
sampling the phase diagram along this one-dimensional 
axis may not reveal the full range of physically relevant 
behavior. A second tuning parameter will be clearly valu- 
able; here we argue that, under suitable conditions, this 
is provided by a magnetic field applied perpendicular to 
the two-dimensional layers. 

Consider the case where both phases in Fig [l] are super- 
conducting; the phase with g < g c then has co-existence 
of long-range order in superconductivity and a secondary 
order parameter. We also restrict attention to the case 



9 The microscopic co-existence of magnetic order and supercon- 
ductivity is not universally accepted but strong arguments in 



its favor have been made recently by Khaykovich et al (2002) 
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where the transition at g — g c is second order (related 
results apply also to first order transitions, but we do 
not discuss them here). Imposing a magnetic field, H, 
on these states will induce an inhomogeneous state, con- 
sisting of a lattice of vortices surrounded by halos of su- 
pcrflow (we assume here that H > H c i , the lower critical 
field for flux penetration). In principle, we now need to 
study the secondary order parameter in this inhomoge- 
nous background, which can be a p roblem of som e com - 
plexity. However, it was argued by Demler et al. (2001) 
and Zhang et al.\ ( 2002| ) that the problem simplifies con- 
siderably near the phase boundary at g — g c . Because 
of the diverging correlation length associated with the 
secondary order parameter, we need only look at the 
spatially-averaged energy associated with the relevant or- 
der parameters. 10 We know from the standard theory 
of th e vortex lattice in a BCS superconductor ( Parks 



whe :e H C 2 3> H c \ is the upper critical field above which 



superconductivity is destroyed. Let us make the simple 
assumption that this change in energy of the supercon- 
ducting order can also be characterized by a change in 
the coupling constant g. We can therefore introduce an 
effective coupling g c g(H) given by 



(12) 



where C is some constant of order unity. As g is linearly 
related to 8, we can also rewrite this expression in terms 
of an effective doping concentration 8 c ff(H), 



(13) 



where C is some other constant. These expressions imply 
that we tune through different values of g or S simply by 
varying the applied magnetic field. 

The sign of C is of some physical importance, and can 
be deduced by a simple argument. It is observed that in 
the lightly doped cuprates, decreasing 8 leads to a sta- 
bilization of an order associated with the Mott insulator 



at t 



lie expense of bite supeicuiidueLiiig order: (Theie is a 
no nmonotonic dependence on 8 from commcnsurability 



effects near 8 = 1/8, but here too the magnetic order is 
stabilized at the expense of superconductivity). As in- 
creasing H clearly suppresses the superconducting order, 
it must be the case that 8 c g(H) decreases with increasing 
H. This implies that C > 0, and indicates a competition 



(Tranquada et al., 1997) between the two ground states, 



or order p arameters, on either side of the quantum criti- 
cal p oint ( Chubukov et al. , 1994a ; gachdev , 2000 ; Zhang 
199%. 



10 Evidence that the primary effect of a magnetic field is a spa- 
tially uniform modification of the magnetic order has appeared 
in recent muon spin resonance experiments 
fjemura, 2003). 
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FIG. 12 Phase diagram in the g,H plane deduced from (p"2|). 
The phase boundary is determined by setting ggg(H) = g c , 
which leads to a phase boundary at a critical field H ~ 
(g — g c )/\n(l/(g — gc))- We assume that g is a monotoni- 
cally increasing function of 8. The collinear magnetic order of 
Figs ^3 and c is the secondary order parameter investigated in 
recent neutron scattering e x perim ents in doped L^CuCU: the 



1961j ) that the energy density ot the superconducting e t a l. 
ord er increases by trie traction ~ (H/H C 2) ln(n C 2/H ), B. T 



observations of |Lake et al. 
those of 



et al 



( 2001 ) arc along the a rrow A, and 
jKatano et al] (|200C ~, [Lake et all ( |2002l ), |Khaykovicb 
( 2002| ), and Khavkovicl i et at. J 2003[) a re a one t he arrow 



l e STM e xperiments of jHoffmanef al\ (|2002 



Hoffman 

fl2002b| ), [Howald et ai\ ( |2002| ), |Howald et qff (p003t > are 



along arrow C, and will be discussed in Section [V.D 



The relationships (|l|) and ([l3]) can be combined with 
Fig to produce a phase diagram in the (g, H) (or (8, H)) 
plane. This is shown in Fig O. Notice that the phase 
boundary comes into the g — g c , H = point with van- 
ishing slope. This implies that a relatively small field is 
needed in the g > g c region to tune a BCS supercon- 
ductor across a quantum phase transition into a state 
with long-range correlations in the secondary order pa- 
rameter. There are also some interesting modifications 
to Fig [l2] in the fully three-dimensional model which ac- 
counts for the couplin g between adjacent C u02 layers; 



these are discussed by Kivelson et al. ( 2002b ) 



A number of neutron scattering studies of the physics 
of Fig [12] in doped La2Cu04 have recently appeared. 
The secondary order parameter here is the collinear mag- 
netic order of Figs ||b and c, which is also observed in 
La 2 - a _5Nd y Sr504 as discu ssed above. Earl i er, a s eries 
of be a utiful expe riments by Wakimoto et al. ( 1999| ), Lee 
et al\ ([L999D , and |Wakimoto et al.\ ( |2001| ) established that 



L^-^SraCuG^ has long-range, collinear magnetic order 
co-existing with superconductivity for a range of 8 val- 
ues above 8 = 0.055. Moreover, the anomalous frequency 
and temperature dependence of the dynamic spin struc- 
ture factor (Chubukov et al, 1994a: Bachdcv and Y(j. 
1992 ) in neutron scattering experiments by Aeppli et al' 



( 1997 ) gave strong indications of a second-order quantum 
phase transition near 8 ~ 0.14 at which the magnetic or- 
der vanished. We identify this transition with the point 
g = g c , H = in Fig |12[ Recent studi es have explored 
the region with H > 0: Lake et al. ( 2001] ) observed a soft- 
ening of a collective spin excitation mode at 8 = 0.163 in 
the presence of an applied magnetic field. We interpret 
this as a consequence of the low H approach to the phase 
boundary in Fig [l2| in the g > g c region, as indicated by 
the arrow labeled A. Notice that the field was not large 
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enough to cross the phase boundary. 

A separate set of experiments have examined the H 
dependence of the static magnetic momen t in the super 



cond u ctor with g < g c in La2-aSr5Cu04 ( Katano et 
200(3; take et all |2002| ) and La 2 CuQ 4+ ^ flKhaykovicb 
eTjl gOO^ " |2002| ) , along the arro w indicated by B in 



2001 



Fig |12j. The theor etical prediction (Demler et al. 
Zhang et al\ , 2002 ) for these experiments is a simple con- 
sequence of (12|) and (|TJ). Let I(H,8) be the observed 



intensity of the static magnetic moment associated with 
the order in Figs ^p,c at a field H and doping 5. If we 
assume that the dominant effect of the field can be ab- 
sorbed by replacing 6 by the effective 5 c s(H), we can 
write 

I(H, 5) « I(H = Q,S cS (H)) 

»/(*-0.fl +*(*;.) !■(**), (14) 



where in the second expression we have used ( |13[) and 
expanded in powers of the second argument of I. Rea- 
soning as in the text below ([l3|) for C, we use the ex- 
perimental fact that a decrease in S leads to an increase 
in the magnetic order, and hence T> > 0. The results of 



recent exper iments ( Khaykovich et al. , 2003 , 2002 ; Lake 
et afy |2002 ) are in good agreement with the prediction 
(14), with a reasonable value for T> obtained by fitting 
(14) to the experimental data. 



B. Detecting topological order 

The magnetic and bond orders break simple symme- 
tries of the Hamiltonian, and, at least in principle, these 
can be detected by measurement of the appropriate two- 
point correlation function in a scattering experiment. 
The topological order of Sections III.B.2 and III.C.2 is 



a far more subtle characterization of the electron wave- 
function, and can only be observed indirectly through 
its consequences for the low energy excitations. We re- 
view here the rationale behind some recent experimental 
searches (Bonn et al, 2001; Wynn et al, 2001) for topo- 



logical order. 

T he searches relied on a peculiar property of a super- 



con ductor proximate to a Mott insulator with topological 
order: there is a fundamental distinction in the internal 
structure of vortices in the superconducting order, spec- 
ified by (|3j> , which depends on whethe r the integer n v is 
even o r odd. This diff erence was noted (Nagaosa and Lcc , 
19921 ; fcachdcvj |l992|) in the context of a simple mean- 



field theory of a superconductor near an insulating spin 
gap state. However, the significance and interpretation 
of the mean-field result, and in particular its connection 
with topological order, did not become apparent until the 
far-reaching work of |Senthil and Fisher] (|2000D, [5cnthil| 
and |Fishcr| ( ^001a| ), and jScnthil and Fisher] ( |2001b| ). The 
arguments behind the dependence on the parity of n v 
are subtle, and only an outline will be sketched here — 



3cnth.il and Fisher (2001b) for a complete exposition. Al- 
though the superconducting order of BCS theory in (||) 
and the topological order of the Mott insulator are quite 
distinct entities, there is an important connection be- 
tween them in the superconducting state: each vortex 
with n v odd in (^) has a vison attached to it. The vi- 
son gap in the proximate Mott insulator then increases 
the energy required to create n v odd vortices, while this 
extra energy is not required for n v even. 

The connection between n v odd vortices and visons is 
most transparent for the case where the spinons in the 
Mott insulator obey fermioni c statis tics. We considered 
bosonic spinons z a in Section III.C.2, but they can trans- 



mute into fcrmions by binding with a vison ( Demler et 



200J; [Kivclson[ [L989fc [Read and Chakrabortyt |1989fl : we 
represent the fermionic spinon by fj a . In the doped Mott 
insulator, each electron annihilation operator, cj a , must 
create at least one neutra l 5* = 1/2 spinon ex citation, 
along with a charge e hole (Kivelson et al., 1987), and we 
can represent this schematically by the operator relation 



(15) 



where &] creates a bosonic spinless hole. In this picture 
of the doped Mott insulator, the presence of supercon- 
ductivity as in (Q) requires both the condensation of the 
bj, along with the condensation of "Cooper pairs" of the 
spinons fj a . We can deduce this relationship from (||) 
and (US) which imply, schematically 



(16) 



where we have ignored spatial dependence associated 
with the internal wavefunction of the Cooper pair (hence 
there are no site subscripts j in ([l6|)), and Af ~ (fjtfj'l) 
is the spinon pairing amplitude. From ( |l6| ) we see if the 
phase of bj winds by 2tt upon encircling some defect site, 
then phase of Ao will wind by A-k, and this corresponds 
to a vortex in the superconducting order with n v = 2 
in (||). Indeed, the only way ( |l6| ) can lead to an ele- 
mentary vortex with n v = 1 is if the phase of the spinon 
pair amplitude, Af, winds by 2tt upon encircling t he vor- 
tex: the latter is another description of a vison ( Benthil 



and Fisher, 2000). This argument is easily extended to 



show that every odd n v vortex must be associated with 
at least an elementary vortex in the phase of Af, thus 
establishing our claimed connection. 

Sufficiently close to the Mott insulator, and near a 
second-order superconductor-insulator transition, the en- 
ergy required to create a vison raises the energy of n v — 1 
vortices, and the lowest energy vortex lattice state in an 
applied magnetic field turns out to have vortices with flux 
hc/e, which is twice the elementary flux ( |5achdev , 1992). 
This should be easily detectable, but such searches have 



not been successful so far (Wynn et all 2001) 



More recently Benthil and Fisher (2001a) have pro 



the reader is referred to Senthil and Fisher (2001a) and 



posed an ingenious test for the presence of visons, also 
relying on the binding of a vison to a vortex with flue 
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hc/(2e). Begin with a superconductor in a toroidal geom- 
etry with flux hc/(2e) penetrating the hole of the torus. 
By the arguments above, a vison is also trapped in the 
hole of the torus. Now by changing either the tempera- 
ture or the doping level of the superconductor, drive it 
into a normal state. This will allow the magnetic flux 
to escape, but the topological order in the bulk will con- 
tinue to trap the vison. Finally, return the system back 
to its superconducting state, and, quite remarkably, the 
vison will cause the magnetic flux to reappear. An exper- 
imental test for this "flux mem ory effect" has also been 
undertaken (Bonn et at, 2001), but no such effect has 



yet been found. 

So despite some innovative and valuable experimental 
tests, no topological order has been detected so far in the 
cuprate superconductors. 



C. Non magnetic impurities 



We noted in Section III.B.l that one of the key con- 



sequences of the confinement of spinons in the bond- 
ordered paramagnet was that each non-magnetic impu- 
rity would bind a free 5=1/2 moment. In con trast, in 
the topologically ordered RVB states of Section III.B.2, 
such a moment is not generically expected, and it is more 
likely that the "liquid" of valence bonds would readjust 
itself to screen away the offending impurity without re- 
leasing any free spins. 

Moving to the doped Mott insulator, we then expect 
no free 5 = 1/2 moment for the topologically ordered 
case. The remaining discussion here is for the confining 
case; in this situation the S = 1/2 moment may well sur- 
vive over a finite range of doping, beyond that required 
for the onset of superconductivity. Eventually, at large 
enough hole concentrations, the low energy fermionic ex- 
citations in the d-wave superconductor will screen the 
moment (by the Kondo effect) at the lowest tempera- 
tures. However, unlike the case of a Fermi liquid, the lin- 
early vanishing density of fermionic states at the Fermi 
level implies that the Kondo temperat ure can be strictly 
zero for a finite range of parameters (Gonzalez-Buxton 
andl lngcrsciit], |l998| ; |Vojta and Bulla| , [2002| ; |Withoff and 
|1990| ) 



Frackin, 1990). So we expect each non-magnetic impu- 
rity to create a free 5=1/2 moment that survives down 
to T = for a finite range of doping in a d-wave su- 
perconductor proximate to a confining Mott insulator. 
The collincar magnetic or bond order in the latter in- 



perconducting critical temperature, and also noted the 
implication of their observations for the confinement of 
spinons, in the spirit of our discussion above. Subsequent 
specific heat a n d nuc l ear magnetic re s onan c e experiments 



(lAlloul et all |1991|; |Bobroff et all |200l| ; |Julicn et 



2000T Sisson et al.[ |200C ) have also explored low temper- 



atures in the superconducting state, and find evidence 
of spin moments, which are eventually quenched by the 
Kondo effect in the large doping regime. Especially no- 
table is the recent nuclear magnetic resonance evidence 
(Bobroff et at, 2001) for a transition from a T — free 



moment state at low doping, to a Kondo quenched state 
at high doping. 

We interpret these results as strong evidence for the 
presence of an S = 1/2 moment near non-magnetic im- 
purities in the lightly do ped cuprates. We have al so ar- 
gued here, and elsewhere ( Sachdcv and Vojta , 2000 ), that 
the physics of this moment formation is most naturally 
understood in terms of the physics of a proximate Mott 
insulator with spinon confinement. 

The creation of a free magnetic moment (with a lo- 
cal magnetic susceptibility which diverges as ~ 1/T as 
T — > 0) near a single impurity implies that the cuprate 
superconductors are exceptionally sensitive to disorder. 
Other defects, such as vacancies, dislocations, and grain 
boundaries, which are invariably present even in the best 
crystals, should also have similar strong effects. We spec- 
ulate that it is this tendency to produce free moments 
(and local spin order which will be induced in their vicin- 
ity) which is responsible for the frequent recent observa- 
t ion of magnetic moments in the ligh tly doped cuprates 



flSidis et al\ |200l| ; [Bonier et al[ [200l| ) 



D. STM studies of the vortex lattice 



Section IV. A discussed the tuning of collincar mag- 



netic order by means of an applied magnetic field, and 
its detection in neutron scattering experiments in doped 
La2Cu04. This naturally raises the question of whether 
it m ay also be possible to detect the bond order of Sec- 



tion 



III.B.l somewhere in the phase diagram of Fig 



Clearly the state with co-existing collinear magnetic and 
superconducting order (explored by experiments along 



sula;or may also survive into the superconducting state, Sachdev, 2001 



but there is no fundamental reason for the disappearance 
of these long-range orders (bulk quantum phase transi- 
tions) to coincide with the zero temperature quenching 
of the moment (an impurity quantum phase transition) . 

A very large number of experimental studies of non- 
magnetic Zn and Li impurities have been carried out. 
Early on, in electron param agnetic resonance experi- 
ments Finkclstcin et al. (199C) observed the trapping of 



the arrow B) should, by the arguments of Section III.C.l , 
also have co-existing bond order. However, more inter- 
esting is the possibility that the BCS supercon ductor it- 
self has local re gions of bond order for H ^ ( Park and 
As we have argued, increasing H in- 



an S = 1/2 moment near a Zn impurity above the su- 



creases the weight of the Mott insulator order parameter 
correlations in the superconducting ground state. The 
appearance of static magnetic order requires breaking of 
spin rotation invariance (in the plane perpendicular to 
the applied field), and this cannot happen until there is 
a bulk phase transition indicated by the phase boundary 
in Fig [l2| In contrast, bond order only breaks transla- 
tional symmetry, but this is already broken by the vor- 
tex lattice induced by a non-zero H. The small vortex 
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cores can pin the translational degree of freedom of the 
bond order, and a halo of static bond order should appear 
aro und each vortex core (Demlcr et al\, |2001 1_ Park and 



Sacl tdcv 



2001 



Polkovnikov et al , 2002a Jq; Zhang et 



2002). Notice that this bond order has appeared in the 
state which has only superconducting order at H = 0, 
and so should be visible along the arrow l abelled C in 
Fig |lj. Recall also our discussion in Section |III.B.1 that 
site charge order is a special case of bond order (with 
r„. = in the bond order parameter Q a (f}). 



Many other proposals have also been made for 
additional order parameters within the vortex core. 
The ea rliest of these involved dynam i c ant iferromag- 
netism( Nagaosa and Lcc , 1992; Sachdev, 1992), and were 
discussed in S ection [1V.B in the context of topological or- 



der. Other s ( Andersen et al. ,_ 2003|; Arovas et al.. 
Chen et al\, |2002|; |Chen and Tirigj, |2002t |Franz et alj 



1997 



2002 



Ghosal et alj |2002fc |Ichioka and Machida] , E002fc |Zhang 



1997; p^hu et al. , 2002) involve static magnetism within 
each vortex core in the superconductor. 11 This appears 
unlikely from the perspective of the physics of Fig 12, in 
which static magnetism only appears after there is a co- 
operative bulk transition to long-range magnetic order, 
in the region above the phase boundary; below the phase 
boundary there are no static "spins i n vortices," but there 
is bond order as d i scusse d above ( Park and Sachdev , 
2001; Zhang et al., 2002). (Static spins do appear in 



the three space dimensional model with spin anisotropy 



and int er-planar couplings considered in Kivelson et al. 
( |2002bD .) A separate proposal involving staggered cur- 



rent loops in the vortex core (Kishine et al., 2001; Lee 



and 



Sha , 2003 ; Lee and Wen . 2001 ) has also been made 

Nanoscale studies looking for signals of bond order 
along the arrow C in Fig [L2^ would clearly be helpful. 
Scanning tunnelling microscopy (STM) is the ideal tool, 
but requires atomically clean surfaces of the cuprate crys- 
tal. The detection of collinear magnetic order in doped 
La2Cu04 makes such materials ideal candidates for bond 
order, but they have not been amenable to STM stud- 
ies so far. Crystals of F^S^CaC^Os+a have been the 
focus of numerous STM studies, but there is little in- 
dication of magnetic order in neutron scattering stud- 
ies of this superconductor. Nevertheless, by the reason- 
ing in Fig [l^, and using the reasonable hypothesis that 
a common picture of competing superconducting, bond, 
and collinear magnetic order applies to all the cuprates, 
it is plausible that static bond order should appear in 
Bi 2 Sr2CaCu208+(5 for large enough H along the arrow C 
in Fig|| 

A number of atomic resolution STM studies of 
Bi2 Sr 2 CaCu208+ < 5 s urfaces have appeared rece n tly (|Hoff- 



et al, 2002a b; Howald et al, 2003, 2002). Hoffman 



et a!j~(2002a) observed a clear signal of modulations in 



[chioka et al 



Also, 

with pre-existing long-range collinear spin order 



(2001) studied the vortex lattice in a state 



the local density of electronic states, with a period of 4 
lattice spacings, in a halo around each vortex core. There 
was no corresponding modulation in the surface topogra- 
phy, implying there is little modulation in the charge den- 
sity. However, a bond order modulation, such as those 
in Figs ^3 and d, could naturally lead to the required 
mod ulation in the l o cal de n sity of states. Othe r stud- 
ies flHoffman et~al[ [2002b| ; |Howald et ajj , |2003| , |2002j ) 
have focused on the H — region: here the modu- 
lations appear to have significant contributions (Byers 
1993| ; Wang and Lcc, 2003) from scattering of the 



et al 



fermionic S = 1/2 quasiparticlcs of the superconductor 
(Sect ion 



, but there are also signals (Howald et 



2003, 20021) of a weak residual periodic modulation in 



the density of states, similar to those found at H ^ . 
Theoretically ( | Howald et~. , 2003; Kivelson et al, 2002a; 
Polkovnikov et al\ , |2003| )7it is quite natural that these 
quasiparticle and or der p arameter modu l ation s co-exist 



Howald et 



fl2002p and |Howald et al\ fl2003|) also pre- 
sented results for the energy dependence of this periodic 
modulation, and these appear to be best modelled by 
modul ations in microscopic bond, r ather than si t e, va ri- 
ables ( |Podolsky etal\ [2003| ; |Vojta| , |2002| ; |Zhangj , |2002j ). 

This is a rapidly evolving field of investigation, and 
future experiments should help settle the interpretation 
of the density of states modulations both at H = and 
H 7^ 0. It should be noted that because translational 
symmetry is broken by the vortices or the pinning cen- 
ters, there is no fundamental symmetry distinction be- 
tween the quasiparticle and the pinncd-fluctuating-order 
contributions; nevertheless, their separate spectral and 
spatial features should allow us to distinguish them. 



V. A PHASE DIAGRAM WITH COLLINEAR SPINS, 
BOND ORDER, AND SUPERCONDUCTIVITY 

We have already discussed two experimental possibil- 
ities for the coupling g in Fig |l|, which we used to tune 
the ground state of the doped Mott insulator between 
various distinct phases: the doping concentration, 6, and 
the strength of a magnetic field, H, applied perpendicu- 
lar to the layers. A simple phase diagram in the small H 
region as a function of these parameters was presented in 
Fig [12], and its implications w ere co mpa red w ith a num- 
ber of experiments in Sections IV. A and IV. D. However, 



even though it is experimentally accessible, the field H 
induces a large scale spatial modulation associated with 
the vortex lattice, and is consequently an inconvenient 
choice for microscopic theoretical calculations. Here we 
follow the strategy of introducing a third theoretical axis, 
which we denote schematically by g, to obtain a global 
view of the phase diagram. As we argue below, infor- 
mation on the phases present as a function of g sheds 
considerable light on the physics as a function of H . 

The crucial role of order parameters characterizing 
Mott insulators in our discussion suggests that we should 
work with a coupling, g, which allows exploration of 
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different ground states of Mott insulators already at 
5 = 0. The range of this coupling should obviously in- 
clude regimes where the Mott insulator has the magneti- 
cally ordered ground state of Fig ||a, found in La2Cu04. 
Now imagine adding further neighbor couplings in (^) 
which frustrate this magnetic order, and eventually lead 
to a phase transitio n to a paramagnetic st ate. 12 As 
disc ussed in Section III.C.l , it was argued ( Read and 
Sac frdevj |l989b| , |199C| ; |Sachdev and Park] , |2002| ) that airy 



paramagnetic state so obtained should have bond order, 
most likely in the patterns in Figs |6|a and b. 

It would clearly be useful to have numerical studies 
which tune a coupling g acting in the manner described 
above. Large scale computer studies of this type have 
only appeared recently. The first results on a quantum 
antiferromagnet which has a spin of S = 1/2 per unit cell, 
whose Hamiltonian maintains full square lattice symme- 
try, and in which it is possible to tune a coupling to 
destroy the collinear magnetic order, were obtained re- 
cently by [Bandvik et al. 




(Gm+B) — ► 



(2002). Their model extended 



with a plaquette ring-exchange term, and had only a 
U(l) spin rotation symmetry. Theoretical extensions to 



this case have also been discussed ( Lanncrt et al 
Park and Sachdev 



2001 



2002 ). Along with the collinear mag- 



nctic state in t he sm all ring-exchange region (small g), 
(Sandvik et al. ( |2002|) found the bond-ordered paramag- 
netic state of Figpk in the large ring-exchange region 
(large g). 

A second large scale computer study of the destruction 
of collinear magnetic order on a mo del with S — 1 /2 p er 



doping 8 



FIG. 13 Zero temperature, zero magnetic field phase diagram 
as a function of the doping 8, and a coupling constant g. 
Here g is, in principle, any coupling which can destroy the 
collinear magnetic order at (tt, 7r) in the undoped insulator, 
while the Hamiltonian maintains full square lattice symme- 
try with spin S = 1/2 per unit cell. The states are labeled 
by the orders which exhibit long-range correlations: collinear 
magnetic (CM), bond (B) and d- wave- like superconductivity 
(SC). At 8 = 0, the CM order is as in Fig |a, the B order is 
as in Fig K)k or b , and we have assumed a co-e xisting CM+B 



unit cell was performed recently by Harada et al. ( [2003 ) . 
They generalized the spin symmetry group from SU(2) 
to SU(iV); in our language, they used the value of N as 
an effective g. They also found the bond order of Fig ^a 
in the paramagnetic region. 

These theoretical studies give us confidence in the the- 
oretical phase diagram as a function nf g and 6. sketched 



r egion , following sachdev and Park (2002) and Sushkov et al 
( |200l|) . The ground state will remain an insulator for a small 
range of 5 > (induced by the long-range Coulomb interac- 
tions), and this is represented by the shaded region. Fhe CM 
order for 8 > could be as in Fig ^b or c, and the B order as 
in Figs b, c , or d , but a variety of other pe riods are also 



possible ([VojtaJ, |2002| [Vojta and Sachdev| , 1999, ). The dashed 
line a indicates the path followed in Fig |12| at H — 0, but 
the physical situation could also lie along line (3. A number 
of other complex phases are possible in the vicinity of the 
m ulticritical po i nt M ; these are not s hown but a re discussed 



in fchang et all ( gOOg ) 



Zaanen et al. 



in dicr p lpHry an H TWl| [Tqqt]- jv^l |9n7^ : ^ Nussinovj (pQOap and also, briefly, in Section |V§ 



( |200l| ), and | aanen and 



and Sachdev, 1999; Vojta et al., 2000a). Phase diagrams 
with related physical ingredients, but with significant dif- 
ferences, ajjjjear in the work of Kivelson et al. (1998) and 
Zaanen| ( [19991 ). 

Important input in sketching Fig [l3] was provided by 
theoretical studies of the effects of doping the bond- 
ordered paramagnetic Mott insulator at large g. In this 
region without magnetic order, it was argued that a sys- 
tematic and controlled study of the doped system was 
provided by a generalization of the SU(2) spin symme- 
try 13 to Sp(2Af), followed by an expansion in 1/N. This 



approach directly gives (Sachdev and Read, 1991) a sta- 
ble bond-ordered state at S = 0, a stable <i-wave super- 
conductor at large 6, and a region in which these two or- 
ders co-exist at small values of 5; all of these phases are 
nicely in accord with the overall philosophy of the present 
article. This analysis of a model with purely short-range 
interactions also found a phase separatio n instability at 
small values of 5 ( Sachdev and Read , 1991 ) , whose impor 



tance had been emp hasized by others ( Bang et al. , 1991 



Emery et al., 1990) on different grounds. With long- 



12 We assume that there is no intermediate state with non-collinear 
magnetic order, as this is not supported by observations so far. 

13 The group SU(2) is identical to the symplectic group Sp(2), but 
the group SU(2iV) is distinct from Sp(2iV) for N > 1. Conse- 
quently, distinct 1/N expansions are generated by models with 
SU(2JV) or Sp(2A0 symmetry. The Sp(27V) choice better cap- 
tures the physics discussed in this article, for reasons explained 



range Coulomb interactions no macroscopic phase sepa- 
ration is possible, and we have to deal with the p hysics 
of frustrated phase separation ( Emery et al. , 199C ). The 
interplay between bond order and d-wave superconduc- 



in Sachdev and Read (1991) 
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ti vity has been studied in some det a il in this reg i on ( Vo- unity of the B+SC phase at H = 



jta |, gOOg ; |Vojta and Sachdcvj , |1999| ; |Vojta et al\ , |2000a|) 



more complex bond ordered structures with large peri- 
ods can appear, usually co-existing with superconductiv- 
ity (as sketched in Fig |l3|). Predictions were made for 
the evolution of the ordering wavevector with <5, and the 
period 4 structures in Figs |^c and d were found to be 
especially stable over a wide regime of doping and pa- 
rameter space. 

The phase diagram of Fig [l3| also includes a region 
at small g, with collinear magnetic order, which is not 



( 


Machida , 


1989 


; Poilblanc and Rice, 


1989; 


1989 


; Zaanen and Gunnarsson 


, 19891)— the 


accu- 



mulation of holes on sites which are anti-phase domain 
walls between Neel ordered regions — is associated with 
this region. However, these stripe analyses treat the mag- 
netic order in a static, classical manner, and this misses 
the physics of valence bond formation that has been em- 
phasized in our discussion here. A related feature is that 
their domain walls are fully populated with holes and 
are insulating. Upon including quantum fluctuations ac- 
counting for valence bonds, it app ears likely to us that 
the stripes will have partial filling (Kivelson and Emery 



1996; Nayak and Wilczek, 1997), acquire bond order, and 



co-exist with superconductivity, as has been assumed in 
our phase diagram in Fig |l3|. Indeed, as we have empha- 
sized throughout, it may well be that the modulation in 
the site charge density — which is proportional to Q a (r) 
with f a = in ( |l0| ) — is quite small, and most of the 
modulation is for r a ^= 0. 

The reader should now be able to use the perspective 
of the phase diagram in Fig pjlto illuminate our discus- 
sion of experiments in Section|V[ The phase diagram in 
Fig |l2|, used t o analyze neutron scattering expe rimen ts in 
Section IV. A and STM experiments in Section IV. D , has 
its horizontal axis along the line labeled a in Fig|l§; the 
phases that appear in Fig [l2] as a function of increasing 
H should be related to those in Fig [D| as a function of 
increasing g, although the detailed location of the phase 
boundaries is surely different. 14 The absence of t opolo gi- 



cal order in the experiments discussed in Section IV. B, 
seen in Fig [l3| to be related to the absence of states with 
non-collinear spin correlations or topological order. The 
formation of S = 1/2 moments near non-magnetic impu- 
rities is understood by the proximity of confining, bond- 



orda red ph. 



in Fig p -3f: — The possible signals of bond 



ord er in a superconduct or at H = in the STM ob acrva 
tions of Howald et al. ( [Howald et al , 2003 , 2002 ), may 
be related to the B+SC phase along the line (3 i n Fig |l3|; 
similarly, the observations of Hoffman et al. (Hoffman 
et a f], p002a| ) at H + 1 can be interpreted by the prox- 



14 More precisely, generalizing the arguments leading to ( Jl 2| ) and 
, we can state that the system is characterized by an effective 
g which increases linearly with H\n(l/H), and an effective 5 
which decreases linearly with H ln(l/H) 



VI. OUTLOOK 

The main contention of this article is that cuprate 
superconductors are best understood in the context of 
a phase diagram containing states characterized by the 
pairing order of BCS theory, along with orders associ- 
ated with Mott insulators; the evidence so far supports 
the class of Mott insulators with collinear spins and bond 
order. The interplay of these orders permits a rich va- 
riety of distinct phases, and the quantum critical points 
between them offer fertile ground for developing a con- 
trolled theory for intermediate regimes characterized by 
multiple competing orders. This approach has been used 
to analyze and predict the results of a number of recent 
neutron scattering, fluxoid detection, NMR, and STM ex- 
periments, as we have discussed in Sections IV. A, IV. B, 
[IV.C , and IV. D. Further experimental tests have also 
been proposed, and there are bright prospects for a more 
detailed, and ultimately quantitative, confrontation be- 
tween theory and experiment. 

All of the experimental comparisons here have been 
restricted to very low temperatures. The theory of 
crossovers near quantum critical points also implies in- 
teresti ng anomalous dynamic propertie s at fin ite temper- 
ature ( Sachdev , 1999 ; Sachdcv and Ye , 1992] ), but these 
have not b een discussed. However, we did note in Sec- 
tion IV. A that the transition involving loss of magnetic 
order in a background of superconductivity was a natu- 
ral candidate for explaining the singular temperature and 
frequency dependence observed in the neutron scattering 
at 5 w 0.14. QAeppli et all |l997j ) 

There have also been several recent experimental pro- 
posals for a quantum critical point in the cuprates at 
S ~ 0.19, linked to anomalous quasiparticle damping 
(|Valla et al , 1999), t hermodynamic ( Tallon and Lora: 



2001), or magnetic (Panagopoulos et al, 2003| , 200 



properties. The study of Panagopoulous and collabo- 
rators presents evidence for a spin glass state below the 
critical doping, and this is expected in the presence of 
disorder at dopings lower than that of the point M in 
Fig H 

Among theoreti cal proposals, a candidate for a quan- 
tum critical point (Sachdev and Morinari, 2002; Zaanen 
et al , 200 1|; Zhang et al, 2002) at large dopings is a novel 



topological transition which can occur even in systems 
with collinear spin correlations. While the topological 
order present in systems with non-collinear spin corre- 
lation leads to fractionalization of the electron (as dis- 
cussed in Section 1V.B), the collinear spin case exhibits 



a very different and much less disruptive transition in 
which the electron retains its integrity, but the spin and 
charge collective modes fractionalize into independent en- 
tities. Note that this fractionalization transition was not 
explicitly shown in Fig [l^, and is associated with an ad- 
ditional intermediate state which may appear near the 



19 



point M. Other theoretical proposals for quant um criti- 
cal po ints are linked to the bon d/charge order ( Kivelson 
et a[} |l998t |Seibold et al\ |1998[ ) in Fi g [y| to order asso- 



ciated with c i rculat ing current loops ( Chakravarty et 
2001] ; |Varma > 1997) which has not been discussed in this 



paper, and to a time-reversal symmetry breaking tran 



sition(Khvcshchcnko and Paaskc , 2001 ; |Laughlin , 1998 



3angiovanni et al, 2001; Vojta et al 
d 



jLa 



OOObl ) between 
id xy superconductors. This last 



- v 2 and d 

x 2 — y 2 

proposal offers a possible explanation of the quasipar- 
ticle damping measurements (Valla et al, 1999). Note 



that this transition does not involve any order associated 
with the Mott insulator. Indeed, the d x 2_ y 2 +id xy order 
can be understood entirely within the framework of BCS 
theory, and experimental support for d x 2_ y 2 + i d xv su- 



per conductivity in re cent tunnelling experiments (Dagan 
and Deutschci , 2001 ) appears in the overdoped regime, 
well away from the Mott insulator. 
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